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Abstract 

Inspired by the AdS/CFT correspondence, we study the pseudoscalar mesons of QCD through 
a dual embedding in a strongly curved extra dimensional spacetime. This model incorporates 
the consequences of symmetry and has very few free parameters, due to constraints from five- 
dimensions and the operator product expansion of QCD. Using as inputs /jr and the pion, kaon, 
and rho masses, we compute the eta and eta prime masses to be 520 and 867 MeV, respectively. 
Their decay rates into photons are also computed and found to be in good agreement with data. 
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1. INTRODUCTION 



While QCD has been unequivocally established as the theory of strong interactions, the 
resolution of one of its mysteries, the U{1) problem, has remained somewhat unsatisfying. 
The U{1) problem is that the Lagrangian of QCD 

^QCD = ^GJ, + cfpq^ + (1) 

has, in the massless limit, a global chiral U{1)a symmetry, under which qi — »• e^^'^^qi, which 
does not seem to be reflected in the spectrum of light pseudoscalar mesons. The formation 
of quark condensates {(jiqj) ~ A^^j^Sij spontaneously beaks the ^7(3)^ x U{3)ji symme- 
try of massless QCD down to a diagonal f/(3)y, which should result in nine pseudoscalar 
pseudogoldstone bosons. The problem is that, with masses included, chiral perturbation 
theory unambiguously predicts a neutral pseudoscalar meson whose mass is strictly less 
than y/Sm-j^ However, the true hadron spectrum contains only the regular 7r° (140), the 
T] (549), and the r]'(957), so the chiral perturbation theory bound is clearly violated. 

Actually, the U{1) problem is little more subtle; the U{1)a is anomalous, i.e. broken by 
quantum effects. Mathematically, while the QCD Lagrangian is invariant (except for the 
mass terms), the functional measure in the path integral is not, and so a chiral rotation 
results in 

^QCD — ^ -^QCD + d — GfivGtiv (2) 
OTT 

where G^u = ^e^^a/sGa/s- The simplest solution of the U{1) problem is then to say that 
the U{1)a symmetry is not really a symmetry at all, so there should be no corresponding 
pseudogoldstone boson [1,[3|. However, from the QCD side it is hard to see how the new term 
in ([2]) could make any difference. Because it is a total derivative, any Feynman diagram 
involving the anomaly carries a factor of zero for total momentum. Thus, the new term 
does not contribute at any order in perturbation theory, and therefore the solution must be 
non-perturbative. 

If we accept that Ui\}A is not a symmetry, the pseudoscalar sector can be modeled in the 
chiral Lagrangian [l, 0, |^ . In full generality, the chiral Lagrangian has five free parameters, 
corresponding to a scale for the U{1) breaking, and four decay constants characterizing the 

interaction strength between the rj and t]' mesons and the jj!^^ and jj^"^ currents 0,111. The 
number of parameters can be reduced, for example by going to the large limit, but a few 
parameters at least remain. In particular, none of the assumptions allow a first-principles 
calculation of the t]' mass. 

The first convincing resolution of the U{1) problem was given by 't Hooft who ar- 
gued topological instanton configurations of the QCD vector potential can contribute to the 
path integral through the anomaly. The instanton contributions are suppressed by factors 
of exp(— l/^fg), which can be significant only for large Qg, i.e. when QCD enters the non- 
perturbative regime. However, instanton calculations generically have infrared divergences, 
due to integrals over large instanton size, so it is impossible to use them for precise quan- 
titative calculations. Nevertheless, they seem to be the correct quahtative solution. And, 
for example, using QCD sum rules 0], they can be used to get a ballpark estimate of the r/' 
mass (~ 1 GeV). 

Other non-perturbative insights into the U{1) problem have come from the lattice js], 0]. 
Because the rj' is critically sensitive to both quark loops and non-local field configurations, 
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it has been a challenge to simulate. Nevertheless, the lattice has been remarkably successful 
in this case, and a recent estimate [8] puts the rj at 871 ± 46 MeV, which is within 10% 
of the experimental value. This is absolute confirmation that QCD itself solves the t/(l) 
problem. But it is hard to get any qualitative understanding from such a purely numerical 
approach. 

In this paper, we propose that the pseudoscalar mesons can be studied both qualitatively 
and quantitatively with a non-perturbative framework based on the AdS/CFT correspon- 
dence jl3]. This framework has already produced an impressive post-diction of the meson 
spectrum, decay constants, and couplings [ll|, 12, It has also led to some new insights 
into observations about QCD, such as vector meson dominance [12], and the structure of 
tensor mesons [iS^. Thus, it is natural to ask whether it can say anything about the f/(l) 
mystery of QCD. 

The approach to AdS/QCD we take in this paper is completely bottom up. Although the 
AdS/CFT correspondence began strictly as a duality between a four- dimensional conformal 
gauge theory and a 10-dimensional string theory, it is difficult to make any quantitative 
predictions about QCD from the string side. Early work, studying for example, the glueball 
spectrum of a large- theory [l3|, had some success when compared to lattice results; but 
to study real world QCD, with three flavors and massive quarks, we would need much more 
information about the string dual of QCD than is currently known (see, for example (l5|). 
Instead, we assume that whatever the string theory is, it must contain bulk modes dual to the 
various local operators of QCD. To reproduce the conformal behavior of the asymptotically 
free regime of QCD, these modes will propagate on a background close to Anti-deSitter space. 
It turns out this is enough information to reproduce a number of non-trivial quantitative 
predictions about QCD at low energy. 



2. SETUP 

The setup is a five- dimensional space, with background metric 

ds" = ^^{dxl - dz^) (3) 

In pure AdS the warp factor is w{z) = 1, but we will allow for background corrections due to 
deviations from conformality. The extra dimension can be thought of as energy, with small 
z representing high energy. Thus we model the IR, where QCD is strong, by boundary 
conditions at a point Zm ~ 1/Aqcd- We also impose boundary conditions at z = 0, high 
energy, where QCD approaches a trivial conformal fixed point. Thus the gravity background 
is modeling energies between Aqcd and infinity. 

For this 5-D description to be equivalent to QCD, it should reproduce QCD correlation 
functions of external currents. These appear as probes in the UV. For each QCD operator 
which couples the vacuum to these currents, there must be a corresponding field in 5D which 
also couples to the current. More generally, for each operator in QCD, there should be a 
5D field. In our case, the currents of interest are J^^ and J^^^, the right- and left-handed 
f/(3) currents. The corresponding fields are bulk gauge fields and Aj^. The operator qiQj 
which spontaneously breaks U{3) x U{3) U{3)v is represented by bifundamental bulk 
scalars Xij. The fields Xjj have interactions and a potential. However, this potential is 
neither calculable nor relevant to low energy, so we simply parameterize this potential and 
fit to data. 
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To study the U{1) problem, we now introduce a new complex field Y to represent the 
square of the gluon field strength: Y ~ G"^^. We can think of the phase of Y as dual to GG. 
We emphasize that identifying Y is not important for the low energy physics, we only use it 
to manifest a linear representation of the symmetries. Thus our 5D Lagrangian, including 
all terms allowed by symmetry, is 



C = + F^) + Tt{\DX\' + 3|Xp} + ^-\DY\' + '^[Y^f det(X) + h.c.]} 



(4) 



That X gets a 5D mass but Y does not follows form the AdS/CFT map between masses 
and dimensions of operators. With these masses, the solutions to the equations of motion 
for X and Y in pure AdS^ are 

(Xij) = Vij{z) = aijZ^ + mijZ (5) 
(Y) = Ez^ + C (6) 

These must correspond to the vacuum expectation values, {(jigj) ~ cTij ~ Aq,;;^ and (G^^,) ~ 
S ~ ^QCB sources Mg and gg. Thus the 2;-dependence of a field is seen to match the 

scaling dimension of the corresponding operator. In the 3-flavor case, for simplicity, we will 
assume that 0",^ = aSij {i.e., (ss) = (dd) = {uu)) and use only two masses m = |(m„ + m^) 
and iTis. 

To study this theory, we will explore the pseudoscalar excitations around the X and Y 
backgrounds. 

Xi, = (X,,)exp(^r7V'') (7) 
Y = (r)exp(m/v^) (8) 

There are of course scalar excitations as well, but these are harder to study as they are 
sensitive to details of the X and Y effective potentials. For the left and right ?7(3) gauge 
fields, we will only need the axial combination A = A^i — Al. The longitudinal modes 
of A^ mix with the pions, so it is helpful to include them explicitly with the replacement 
d^ip. Then we get 

flavors 

+ ^[(^.^° - d,^f - (^5 - d.af] + ^^^^(« - (9) 

The on r^^ refers to the r° = ;^diag(l, 1, 1) generator of f/(3), and an 8 superscript will 
refer to the = ^^diag(l, 1, —2) generator, in the u, d, s basis. Note that we have absorbed 



a constant into the definition of C, and absorbed factors of C into the definition of k. We 
have also dropped H ~ {Gl^J)) as it is will be a subleading power correction in everything 
that follows. Regardless of these conventions, it is simplest to regard Eq. ([H]), instead of 
Eq. (jlj), as the starting point for phenomenological analysis. 

In 4D the U{1)a symmetry is anomalous; it is broken by quantum effects. But quantum 
effects in 4D correspond to classical effects in 5D, so the symmetry should be explicitly 
broken in 5D. In unitary gauge, this is true. But in the form ([9]), we have restored the 
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symmetry with our "axion" Goldstone boson a. In fact the whole U{?))a is gauged, so there 
is a local symmetry under which 



a a + a° 

We can use this to set = 0, however it is helpful to retain these modes to simplify the 
calculations. 

Because of the gauge symmetry, the fields ip,!] and a are not strictly independent, but 
they do have different physical meanings as can be seen by introducing external sources (and 
notation). We define 

■/J^Ews^^. (10) 



The U{1)a current is normalized as 



/2Nf 

V J I 

A source J^y4^5(z) on the UV brane leads to {p^dfj_J^S{z) after introducing (p and integrating 
by parts. So the source for ip is 

4 ^ (12) 
Finally r]'^ and a by definition correspond to specific 4D fields, so we have 

4 = 9vQt^5rij<lj^ Ja = ga^GG (13) 



Note that we use constants g^^ and ga to normalize and Ja, while the normalization of 
is set by the interaction strength g^, in the Lagrangian. 

These currents help us identify our pseudoscalar fields. We see that although ip, r] and a all 
mix they still have physical meanings: for a particular mode, a is the "glueball" component 
and 1] and ip are the "quark" components of the corresponding mesonic wavefunction, with 
ip related to the longitudinal mode of the axial vector field. 



A. Matching to QCD 



We will now calculate the parameters in our model by matching to the QCD operator 
product expansion (OPE). Let us start immediately with the case of interest, 3-fiavors, 
massive quarks, and a physical t]'. We will need to make use of the anomaly equation 



4 = 8,4 = VWf^GG + 



(14) 



J flavors 



The anomaly shows up in the OPE |16l . Il7 



{44) 



167r4 



QnogQ' + 



(15) 



flavors 
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As explained in detail elsewhere ll|, ll3|] , correlations functions are calculated in 5D by 
solving the equations of motion in the presence of a source. For example, this lets us deduce 
that (75 = 27r in the case of interest, Nq = 3. For the ip^ correlator, which is relevant for the 
anomaly, we can write 

Here, ^^{z) is a bulk-to-boundary propagator, that is, a solution to the equations of motion 
with (/3°(0) = 1. For this calculation, chiral symmetry breaking is irrelevant to leading order 
and so we can set f = 0. The y?" and equations of motion then become 



d)-d,^' - gl'^{^' - a) = (17) 
glC^d^a - Q'^z^d^if' = (18) 



These are solved perturbatively near 2; = by 



/ = 1 - gl^ \og{Q'z') + al'^^'Q' ^og{Q'z') + ■■■ (19) 
a = ~Q^z^ + ■■■ (20) 



Matching ([16]) to ([T5]) leads to 



a. 



^ = ^V2A^/ (21) 

Note that for this matching we have assumed that as is constant in the UV. Of course, 
as runs with scale, and it is therefore reasonable to assume that as would be a function of 
z as the 1-loop QCD jS function. Hence, we should take 

^=v^^' "^ = 7rr^^ V l3o = ^i^Nc-lNj) (22) 

27r^ /5o log(AQCD2:) 27r 3 3 

where Aqcd ~ z^^. The fact that as varies slowly in the UV, and that d^as ~ a^, makes 
the above matching correct to leading order in a^- 

Instead of sourcing ip^ with 9^ J° 7^ 0, we can also consider pure gluodynamics and source 
a by turning on GG. This will fix the normalization of J^. The QCD correlation function 
of interest is [3] 

2 

XtiQ) = ((^GG){^GG)) = -^Q'logg^ (23) 

which should match 

Xt(g) = AaJ.) = ^lim^ (24) 

gi gi z-^ 

For a solution with a(0) = 1. Solving the equations of motion perturbatively 

a = 1 + ^Qh' - ^Q'z^ log Q^z^ + ■■■ (25) 

lets us deduce that 

Qa = 27r^— = yfWf (26) 
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We will use this later on to compute the topological susceptibility Xt{0). 

Next, we would also like to consider modifications to the background due to deviations 
from conformality, in particular the effect of the strange quark mass. Consider the transverse 
part of the axial vector OPE 



(J'-Ju) = iQ.Qu - v,uQ'')nA + ■■■ (27) 



IIa = -— logg2 - ^M,(gg) + • • • (28) 



The logQ^ term on the right hand side is determined by conformal invariance and is used 
to fit (75. The second term is the power correction in which we are interested now. By 
dimensional analysis, this should modify the warp factor to 

w{z) = 1 + C4Z^ C4 ~ Mg{qq) (29) 

This warp factor modifies the axial-vector equation of motion to 

z z z-^ 

We can solve this perturbatively in C4 ~ mga. For rriq = a = 0^ the solution is 
A^{z) = QzKiiQz). Then, using the AdS Green's function K\z^z') [18], the perturba- 
tive inhomogeneous solution can be written as 

A^{z) = / dz'K{z,z')[2glmqaQz^lCi{Qz) +AciQ'^z^lCQ{Qz)\ (31) 
2q'im(,a Aca\ z^ o, ,„ , 



Thus AdS gives 

= lim^-;-::^:^ = ^logg^ - ( + ^ + . . . (33) 



^9^^2_ , Um^G 8c4\ 1 



Comparing to fl28|) . we deduce 5^5 = 27r and 



C4 = "^gliUq^qq) - ^m^a) = -ym^a (34) 



The last equality follows from the Gell-mann-Oaks-Renner relation jll 



B. Fitting to data 

Having determined g5,C,ga and C4 from matching to the OPE, our Lagrangian is com- 
plete. The remaining unknowns must be fit to data. We use 

rup = 770 MeV z~^ = 323 MeV = Aqcd (35) 

U = 93 MeV a= (333 MeV)^ (36) 

= 140 MeV m = 2.22 MeV (37) 

rriK = A9A MeV = 40.0 MeV (38) 
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FIG. 1: Masses of rj and rj' as the parameters k and Aqcd ai"e varied. In the k plot, we fix 
Aqcd = 320 MeV « z"^ and in the Aqcd plot, we take k = oo. Dashed lines are the experimental 
values. 



This gives C4 = —0.676 for strange and C4 = —0.037 for up/down. Thus we are justified in 
only turning on C4 for the strange quark. Keep in mind that although these "quark masses" 
may seem small, care must be taken when comparing them to masses deduced from another 
scheme. 

There is one more parameter in our Lagrangian that remains, k. The k term corresponds 
to an entirely non-perturbative effect. However, it multiples a function which grows like 
SO we expect it to act effectively like a boundary condition forcing a{zm) = V^izm)- 
Thus, we leave k as a free parameter and show that for the rj' the results are fairly indepen- 
dent of K for large k. 



3. THE 77' 



Having fit all the parameters in our Lagrangian (using only nip, /^r, and rriK), we can 
now look at what masses are predicted. For the neutral pseudoscalar spectrum, there is a 
competition between the quark masses, which force the t] and r]' into the q and s bases, and 
the anomaly and k terms, which push towards the 77° and t]^ basis. There are seven fields, 
^0,8^^0,8^^0,8 ^j^g residual gauge invariance to set A^'^ = 0. Thus, to 

determine the t] and rj' masses, we need to solve a set of five coupled differential equations. 
We find the smoothest numerical results if we use the equations of motion for a, cp'^, cp^, A\ 
and Ag, where 




The equations are 

d.'-^d.^^ - gfl{^^ - ,n + 9l^^{^|\^^ - i=^^ - a) = (42) 
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glCd^a = 



(43) 
(44) 



with Vq = mqZ + az^, Vg = + az^ and C given in equation (l22l) . All the modes have 
Dirichlet conditions in the UV and Neumann in the IR. To canonically normalize the fields, 
we demand 



dz 



r^2 



(45) 



The resulting masses are shown as a function of k on the left side of Figure 1. These 
curves are convergent, and the asymptotic values for large k, as compared to the experimental 
central values (in MeV) are 



nir, 



520 
867 



[m: 
(m: 



EXP 

V 

EXP 
V' 



549) 
957) 



(46) 
(47) 



So we are off by 5% and 9% respectively. We can also turn off the anomaly by lowering 
Aqcd, as shown on the right in Figure 1. 

It is worth emphasizing that taking k ^ oo does not send m^/ oo. In the chiral 



Lagrangian, there is a parameter like k which should be proportional to the anomaly 
and provides a mass term for the U{1) pseudoscalar. In that case, taking k —>■ oo does 
decouple the r]', and the correct t]' mass can only be reproduced by tuning k against the 
other chiral symmetry breaking terms in the Lagrangian. In AdS, we could have simply taken 
K = oo to begin with, which would be a simpler model with k, is replaced by a boundary 
condition. However, we choose to allow k to vary because it gives us an additional handle 
on the U{1) sector. 

Next, we calculate the decay constants. There is not a single and /^/. Instead, there 
is a decay constant for each into the J° and J* currents. 



(48) 
(49) 



These can be calculated from the wavefunctions directly, using relations similar to those 
in For example, we solve the above differential equations with m = rrir,, then evaluate 



fr]0 
frjS 



1 

hm 

gi z 

1 d,^' 
lim 

gi z 



17.0 MeV 
103 MeV 



Similarly 



fv'o 

fri'8 



129 MeV 
-35.1 MeV 



(50) 
(51) 



(52) 
(53) 



So qualitatively, the rj' is more 7]^ and the 1] more rj^, as expected. These values can be 
compared to decay constants extracted within chiral perturbation theory [l^. It is, however. 
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FIG. 2: Profiles of the bulk wavefunctions of the components of r] (left) and rj' (right), for n = 20. 
Because of the z-dependence, there is no simple mixing-angle interpretation. 



misleading to represent this mixing in terms of angles because the \ri^) and \ri^) components 
of the mass eigenstates |?7) and {r]') depend on z. This can be seen from Figure [21 which 
shows the profiles of the AdS wavefunctions of r] and rj' . 

These decay constants are not directly observable. What is observable are the neutral 
pseudoscalar decays P — 77, which are mediated by the axial anomaly. Amusingly, the 
form of this anomalous interaction in five dimensions was derived long go by Wess, Zumino 



and Witten (WZW) [20|, l2l|]. The bulk Chern-Simons (CS) term relevant for the decay to 
photons is 

Ccs = ^e^^c°^y^^l^cDA^^Tr[QV^] (54) 

where Q is the generator of electric charge. Here, Vab are components of the field strength 
for the vector gauge field, Vm{z), from which we want to extract the constant photon zero 
mode by setting V^{z) = 1 (this normalization is consistent with ( fTTI) . see IJ] for more 



^wzw = —e^'""'V^,Vparf"Yx\Q\''\ (55) 



details about the photon). In addition to this bulk term, there is a WZW term on the IR 
boundary at z = z™, 

27r2 

which absorbs the anomaly. For constant V^, with = as usual, the CS term is a 
total derivative, and therefore only the boundary WZW term contributes. Explicitly, the 
amplitude is 

^ " 1 . 4 



(56) 



which leads to (as compared to the experimental values extracted from the observed decay 
rates), in units of TeV~^ 

A,,, = 24.3, {A^f^ = 24.9) (57) 
Ar,'^^ = 48.1, (Aj,^P = 31.3) (58) 

These are the asymptotic values at large k. The variation of the decay constants with k, and 
Aqcd is shown in Figure [31 

We can also find the value of k, which provides the best fit to the experimental val- 
ues of mj^,mri', Ar^^^ and Ar^i^^. This is given by k = 26.1 with = 466 MeV, m^/ = 
792 MeV, A^^^ = 30.2 TeV"^ and A^/^^ = 37.3 TeV"^ The RMS error is 18%. 
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FIG. 3: Decay amplitudes for r] and rj' as k and Aqcd are varied. Dashed lines are the experimental 
values. 

4. TOPOLOGICAL SUSCEPTIBILITY, INSTANTONS, AND 



Now, let us turn to the topological susceptibility, Xt- The standard argument is that 
if there are massless quarks in the theory, then 6 is unphysical, and thus Xt must vanish. 
However, if all quarks are massive, or there are no quarks at all, then we expect Xt to be 
nonzero. These facts lead the Witten-Veneziano relation (22I . 23] for the rj' mass at large Nq 



{w? + mi, - 2m|-) 



(59) 



This relation, which gives Xt = (171 MeV)^, is only approximate. It assumes all the decay 
constants are equal, that the mesons have no glueball component, and that Nc is large. Nev- 
ertheless, lattice seems to confirm these approximations by producing xt = (191MeV)^. 
With our 5D construction, we can calculate the topological susceptibility, the meson masses, 
and the decay constants directly, and furthermore we can verify that Xt vanishes only with 
massless quarks, from which the Witten-Veneziano relation follows. 
Recall that 

Xt = — lim — ^ (60) 



9l 



^0 z-^ 



for a solution with a(0) = 1. First, consider the case of pure gluodynamics. Then there are 
no 7] or fields, and the equation of motion at zero momentum is simply 



dz-Kdza = 



(61) 



In the absence of a k term, it is simplest to just impose a{zm) = directly. Then the 
solution is a{z) = 1 — {^)* = 1 — jidl/ C'^)Xtz'^ from (160|) . Now suppose there are quarks. 
In the limit that the anomaly is weak (for example at large Nc), we can do a perturbation 
expansion in C. To leading order in C, the rj — ip system decouples from the a mode. Then 
from fl50l) we get (p{z) = 1 -|- \g1frfZ^. Then the equation of motion ffTSjl . with Q"^ = —m^ 
and using fl2Tl) and fl26|) . gives 



Xt 



AN, 



r2 2 



(62) 



which matches Witten-Veneziano. 
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To see that xt vanishes with massless quarks, we no longer assume that C is small. Then 
the ?7 and equations of motion (at Q = 0) are 

(mqZ + az'^f . . 

dz — - — ^ dzT] = (63) 

z-^ 

{mqz + az^dzV - C^d^a = (64) 

If rriq = then the only solution for rj satisfying ri{0) = is ri{z) = 0. Then a{z) must be 
constant and the topological susceptibility vanishes. However, as long as niq ^ 0, there is a 
solution with t] ^ z"^ and a ~ z'^ near z = 0. In this case Xt is nonzero. In fact, we can solve 
the equations exactly for constant C with r]{zm) = a{zm) boundary conditions, giving 



vr"^ + niaZmimaZm. + crz^ ^ ™ 



4 

- ■/ 

m 

m) 



For C = this reduces to the result from pure gluodynamics. If C 7^ 0, then we can see 
directly that = forces Xt to vanish, as expected. 

Note that we have not used the k term at all to calculate the topological susceptibility; 
we have only used the fact that it leads to a{zm) = in pure gluodynamics, or a{zm) = 
ri{zm) if quarks are included. The k term is supposed to represent some non-perturbative 
effects which are normally associated with instantons, so it is natural to ask if we can make 
the connection more precise. In QCD the one instanton contribution to the topological 
susceptibility can be calculated explicitly 0, 113] 

xm = ----l rdp^[QV/C2(Qp)]' (66) 
Jo P 

Here, V^p) is the dilute-gas instanton density. For example, for Nc = 3, T>{p) = (AqcdP)^^- 
The Bessel function /C2 appears as the Fourier transform of GG evaluated on a one-instanton 
solution. This expression is divergent due to large instantons, so one normally cuts off the 
integral at p = pc ~ ^qcd- 

In QCD it is not meaningful to compare the contribution of this specific gauge configura- 
tion to any particular calculation on the AdS side. This is because the five-dimensional dual 
describes only gauge invariant quantities resulting from integration over all gauge configura- 
tions, and it is not clear in which sense this particular configuration dominates the integral. 
Nevertheless, in truly conformal theories, where the coupling constant is a marginal param- 
eter, it makes sense to compare non-perturbative contributions to correlators {i.e. in powers 
of e~^/^^) between the CFT and the five-dimensional theory. Of course, the axion would 
have a similar bulk description in such a case, the main difference being the 2;-dependence 
of the K-like term in the dual to the CFT. We thus expect that the n term contribution to 
some correlation function to have structures similar to those one gets from integration over 
instanton size, since the bulk integration over the 2;-variable, must ultimately reproduce the 
same correlation function in the CFT. 

To see the similarity to the instanton calculation, let us look at the GG two-point function 
from the bulk perspective. We can solve for a perturbatively around the conformal limit. 
In the conformal approximation, there is no IR brane, and the axion bulk-to-boundary 
propagator (with a(0) = 1) is 

^35^J_9^a(o)_g2^(o)^0 ^ a^'\z) = lz'Q'}C2{Qz) (67) 
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Conformality is broken by the IR brane and by the k term. With these effects, the equation 
of motion becomes 

C^z^d,\d,a - C^Q\ + ^K(z)a = (68) 



z^ 



where k{z) includes the non-conformal 2;-dependence of the k term. Now, we can get a 
simple expression for the topological susceptibility by integrating the action by parts on the 
equations of motion 



dz 



(J2 (J2 I 

-^{dzof H — ^Q'^q^ ^K(z)a'^ 

^.i ^.i 



CMim^ = 2iV,xt(Q) (69) 



If conformal invariance is a good approximation, we can estimate the effect of conformal 
symmetry breaking by evaluating this expression on a*-°^ We thus find 

Xt{Q) = [Q'z'lCiQz)] ' (70) 

This has exactly the same form as the instanton contribution. Thus, the scale dependence 
of the K term acts just like the instanton density and the IR brane provides a natural cutoff 
on the integral over instanton size. 

Finally, let us say a word about the QCD vacuum angle 9. This angle is intimately 
tied to the solution of the U{1) problem. The argument, roughly, is that the topological 
susceptibility must be nonzero to spht the rj' from the t] and the 7r°. Since the topological 
susceptibility is the second variation of the effective action with respect to 6, there must be 
sensitivity to 9 in QCD. Thus the strong CP problem, which is why the apparent value of 6 
is so tiny {6 ^ 10~^), must be taken seriously. 

In AdS, there are three angles, appearing in the X and Y vevs, and in the k term. We 
can write {X) = |(X)|e'^% (Y) = \{Y)\e^^^ and k = |K|e*^^. Since 6i and 62 come from vevs, 
they can be functions of z (as in Eqs.dSED), but 6^3, like k, should be a constant. In full 
generality, 61 can have flavor indices as well. This leads to 

^ = ^[^5 + d.iv' - OD? + ^[Al + d.{a - e^)]' + ^y^ia - r,' - (71) 

Now, the axial symmetries in AdS are local gauge symmetries, so we can rotate 9\ and 62 
into yf and a respectively. This leaves ^3 = ^ as the physical vacuum angle. Although the 
combination a — yf couples directly to 9, it cannot be the physical axion which solves the 
strong CP problem. Even though (a — rf) = 6,6 cannot be eliminated since it is a and rj 
which appear in the rest of the Lagrangian, not the orthogonal combination a + 77°. However, 
there is hope that since AdS allows a quantitative study of confinement, a strong-dynamics 
based solution to the strong CP problem might be realizable. 



5. SUMMARY AND CONCLUSIONS 

We have studied the U{1) problem through and extra-dimensional model inspired by the 
AdS/CFT correspondence. This model is built from the bottom up, by fitting some param- 
eters to perturbative QCD correlation functions and others to data. All of the parameters 
in the model can be determined by the experimental masses of the vr", and p mesons. 
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and the pion decay constant f-^. This is only one more experimental value than is needed to 
define QCD itself (in QCD, we have the quark masses rriq and and the value of Aqcd)- 
In particular, no strong-dynamics based observable, such as the topological susceptibility, 
is needed to study the rj'. Instead, we only need the coefficient of the anomaly which is 
perturbatively calculable and 1-loop finite. The non-perturbative effects are represented in 
our model with a k, term, on which we have shown the observables are only weakly depen- 
dent. Using this construction, we have calculated m^/ = 867 MeV, which is 9% off from 
experiment. We have also calculated its decay constants, and its coupling to photons, as 
well as the analog quantities for the t]. The best fit for k matches the four observables to 
18%. 

In addition to being quantitatively precise, the extra-dimensional construction allows for 
additional qualitative insight into the U(l) problem and related issues. For example, we 
have shown how the vanishing of a quark masses would cause the topologically susceptibil- 
ity to vanish, independent of any discussion of the theta angle 6 of QCD. From this, the 
Witten-Veneziano relations follow. In QCD, it is difficult to study the contribution of non- 
perturbative effects, because one cannot turn off the anomaly except by taking A^^^ — >■ oo. 
In the holographic model there are two additional parameters, k and Aqcd which can be 
separately dialed, giving us new handles on the anomaly. We also showed that the non- 
perturbative contribution to the topologically susceptibility, which can be represented with 
an instanton calculation, has a direct analog in AdS. The same Bessel functions appear in 
both cases, and the integral over instanton size is replaced by an integral over the extra 
dimension. Instead of having to invoke a separate cutoff to regulate the IR divergence, we 
naturally use the same IR cutoff we would have in a non-anomalous theory. 

This solution to the U{1) problem demonstrates the versatility of the bottom-up 
AdS/QCD approach. It also emphasizes that AdS is not just a complicated way of phras- 
ing the predictions of chiral perturbation theory - the rj' mass is simply a free parameter 
in the chiral Lagrangian. Although our effective description is non-renormalizable, higher- 
dimension operators are quantitatively irrelevant for the observables in question, as is ex- 
pected from naive dimensional analysis. It is therefore likely that through further application 
of the AdS/QCD correspondence, additional quantitative and qualitative information about 
the non-perturbative structure of gauge theories can be derived. 
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